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Abstract 

We consider a family of random line tessellations of the Euclidean 
plane introduced in a much more formal context by Hug and Schneider 
[Geom. Fund. Anal. 17, 156 (2007)] and described by a parameter 
a > 1. For a = 1 the zero-cell (that is, the cell containing the origin) 
coincides with the Crofton cell of a Poisson line tessellation, and for 
a = 2 it coincides with the typical Poisson- Voronoi cell. Let p n {ot) 
be the probability for the zero-cell to have n sides. By the meth- 
ods of statistical mechanics we construct the asymptotic expansion of 
log p n (a) up to terms that vanish as n — > oo. In the large-?! limit the 
cell is shown to become circular. The circle is centered at the origin 
when a > 1, but gets delocalized for the Crofton cell, a = 1, which 
is a singular point of the parameter range. The large-n expansion of 
logp n (l) is therefore different from that of the general case and we 
show how to carry it out. As a corollary we obtain the analogous 
expansion for the typical n-sided cell of a Poisson line tessellation. 

Keywords: random line tessellations, Crofton cell, exact results 
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Figure 1: Example of a line tessellation of the plane. The projection of the 
origin onto line a is denoted R a ; this projected point fully determines a. The cell 
surrounding the origin (the 'zero-cell') has been shaded. 

1 Introduction 

1.1 Poisson-Voronoi and Poisson line tessellations 

A Voronoi tessellation of IR 2 is a partitioning of the plane into cells con- 
structed around 'point particles' in such a way that each point of space is 
in the cell of the particle to which it is closest. When the point particle 
distribution is uniformly random, (or, in mathematical terms, when it cor- 
responds to a homogeneous Poisson point process), the resulting partition is 
said to be a Poisson-Voronoi tessellation. It is one of the simplest mathe- 
matical models of naturally occurring planar cellular structures. Because of 
the great variety of their applications, the statistics of Poisson-Voronoi cells 
has been studied in many different areas of science. References may be found 
in the encyclopedic review by Okabe et al. [I]. 

A different way of partitioning the plane into cells is by means of inter- 
secting straight lines, as in figure [TJ When these have a uniform distribution 
(or, in mathematical terms, when they correspond to a homogeneous Poisson 
line process), we refer to the partition as a Poisson line tessellation. The 
Poisson-Voronoi tessellation and the Poisson line tessellation are both statis- 
tically invariant under translations and rotations in the plane. For both, the 
cells are convex polygons. 

An early application of the Poisson line tessellation occurs in work by 
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Goudsmit [2], carried out at the suggestion of Niels Bohr. The question in 
that paper arose from cloud chamber experiments: when three lines seem- 
ingly originate from the same point, then what is the probability that they 
do not result from the same event? Hence the problem became to calculate 
the probability for three independent lines to nearly pass through the same 
point, or, put differently, for a typical triangular cell to have an area less 
than A in the limit of very small A. Poisson line tessellations have since 
interested mathematicians with important contributions due to, in particu- 
lar, Miles PHE], Matheron [6], Kovalenko Ej, Goldman [9], and Hug, 
Reitzner and Schneider [TO] . 

The cell that contains the origin is generally called the 'zero-cell'. In 
the Poisson line tessellation it carries the special name of 'Crofton cell', in 
reference to Crofton's formula in integral geometry [UJ . Since the origin falls 
in a cell of area A with a probability proportional to A, the Crofton cell is 
not typical but more likely larger-than-typical. In mathematical terms, the 
density of the typical cell differs from that of the Crofton cell by a factor 
A/(A) typ , where (A) typ is the average typical-cell area. We will distinguish 
quantities pertaining to the typical cell of a Poisson line tessellation (as 
opposed to the Crofton cell) by an extra superscript 'typ'. 

1.2 Cell sidedness 

The cell property most studied is the sidedness probability p n , that is, the 
probability for the cell to have n sides. Other quantities of interest have 
included the averages, moments, and correlations of n, the cell area and the 
cell perimeter, as well as the distribution of the angles between the perimeter 
segments. The statistical properties of an n-sided cell may be expressed an- 
alytically as 2n-fold integrals on the planar coordinates of the point particles 
(for the Voronoi tessellation) or of the lines (for the line tessellation) defining 
that cell. However, only few of these integrals can be evaluated exactly. In 
particular the calculation of the fraction p n has so far been impossible for 
general n, whether for the Poisson- Voronoi or for the Poisson line tessellation. 

We briefly recall some known results. In a Poisson line tessellation the 
typical cell has an average number of sides (n) typ = 4. The sidedness p typ 
peaks at n = 4. Miles [3] obtained the exact value p^ p = 2—n 2 /Q = 0.35506... 
and Tanner [TJ showed that pf p = 7r 2 log2 - § - ^vr 2 - |C(3) = 0.38146.... 
Numerical values for p typ based on Monte Carlo simulation were given by 
Crain and Miles [13] , by George p3] , and recently by Michel and Paroux [T5] 
(who compare their results to the earlier ones), for sidednesses not exceeding 
n = 12. 

For the Crofton cell, by contrast, Matheron [6] showed that the aver- 
age sidedness is (n) = n 2 /2 = 4.9348.... The distribution p n peaks at 
n = 5. Miles [5] obtained the only known exact result, namely p-^ = (25 — 
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361og2)7r 2 /6 = 0.076820.... Numerical values for p n were given by Calka [16] 
for n = 3,4, ... ,9 and by Michel and Paroux [15] for n = 3, 4, . . . , 11. 

The Poisson-Voronoi cell has an average sidedness of (n) =6. Its sided- 
ness distribution p n peaks at n = 6. Although no exact results are known for 
any of the p n , it has recently been possible [T71 [TBI dH] to obtain the asymp- 
totic expansion of p n in the limit of asymptotically large n. The expansion 
was shown to have implications for the asymptotic cell shape, as well as for 
the finite n behavior [20] and for correlations between neighboring cells [21] . 

A natural question to ask, then, is whether a similar asymptotic analysis 
of p n (for the Crofton cell) and of p^ p (for the typical cell) can be carried out 
for the Poisson line tessellation. Now it was observed by Hug and Schneider 
[22] that the Crofton cell and the typical Poisson-Voronoi cell are particular 
instances of a more general family of zero-cell problems, dependent on a 
parameter a. Hence our interest in the Poisson line tessellation leads us 
quite naturally to study here this full one-parameter family, which we will 
define in the next subsection. 

1.3 A family of tessellations 

Let a be a line in the plane and R a the projection of the origin onto that 
line, as shown in figure [TJ then a is uniquely specified by R a . In the Poisson 
line tessellation the R a are independent identically distributed stochastic 
vectors with a density proportional to 1/R a - This tessellation is statistically 
invariant under translations. In fact, the intersections of the lines a of the 
tessellation with an arbitrary additional line constitute a one-dimensional 
Poisson process; and the associated angles of intersection 6 a are mutually 
independent and have the common probability density \ sin 9, where < 
6 < 7r. In this work we will consider a more general tessellation that was 
introduced by Hug and Schneider [22] and depends on a parameter a. The 
projection vectors of this tessellation are distributed with a density 

p(R) = est x R a ~ 2 , a>l. (1.1) 

For generic a the distribution (11 .ip has central symmetry around the origin 
but the corresponding tessellation is not translationally invariant. The zero- 
cell is therefore unlike any other cell. For a = 1 we recover the Poisson line 
tessellation. For a = 2 the zero-cell corresponding to (11.11) is identical to 
the typical cell of the Poisson-Voronoi tessellation, which is easily seen as 
follows. For a = 2 we have p(R) = est, so that the projections themselves 
have a uniform density. Then, for a given configuration of projections {R a }, 
we may imagine point particles located at the set of positions {2R a }, as well 
as an extra point particle placed in the origin. These particles constitute a 
Poisson point process in the plane. The Voronoi cell of the particle in the 
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origin is then equal in distribution to the typical Poisson-Voronoi cell [23] . 
while also being identical to the zero-cell of the a = 2 line tessellation. 

Although not at present of any known application in physics, it is instruc- 
tive to study the entire a dependent family of problems in order to see how 
it links together the two cases of greatest renown, the Crofton cell for a = 1 
and the typical Poisson-Voronoi cell for a = 2. 



1.4 Results 

For the a dependent tessellation defined by (11.11) we construct the large-n 
expansion of the probability p n (ct) that the zero-cell is n-sided. In section 
121 we consider the parameter range a > 1, for which our work is a rather 
straightforward extension of earlier work [18] on the Voronoi cell, a = 2. 
The final result for p n (a) is best expressed with the aid of the auxiliary 
quantity pn\a) given by 

P^H^gp, a>l. (1.2) 

In section [21 we derive that when n — > oo, 

p n {a)~C(a)p®{a), a > 1, (1.3) 

where the symbol ~ is defined by the classical equivalence (u n ~ v n if and 
only if u n /v n — >1 for n — > oo), and where the prefactor C(a) is given by 




Equivalently, we have that when n — > oo, 
logp n (a)= - 2n \ogn + n log(vr 2 e 2 a) - \ logn + log f ^g^~ J 



+ o(l). (1.5) 



For a = 2 this result reduces to the known expressions [18] of the Voronoi 
problem. The study for a > 1 shows, however, that the Crofton cell corre- 
sponds to a singular point at the lower limit, a = 1, of the parameter range 
studied here. We have lim^i C(a) = oo, the divergence being due to the 
factor of index q = 1 in ( 11. 41) . This signals the breakdown of the method of 
section [21 when applied to the Crofton cell, a = 1. 

In section [3] we develop the modified approach necessary to find the sid- 
edness probability p n (l) of the Crofton cell, as well as its other statistical 
properties. The appropriate result for the Crofton cell derived there is that 
for n — > oo 

p n (l) ~ frrf>(l), (1.6) 
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or equivalently 

logp n (l)= - 2nlogn + nlog(vrV) + |logn - log(6vr 5/2 ) + o(l). (1.7) 

As a corollary we find in section [4] that the typical cell in a Poisson line 
tessellation has a sidedness probability p^ yp (l) given by 

p*P(l) ~ fn-^a), ^-oo, (1.8) 

The ratio p n (l)/»^ yp (l) ~ |n 2 is that of the area of the re-sided Crofton cell 
to the area of the average cell in a Poisson line tessellation. 

The derivation of these asymptotic series for the p n (a) is based on a 
perturbation expansion around the regular n-sided polygon centered at the 
origin. The prefactor C(a) represents the partition function of the elastic 
deformations of this n-gon, the 'elasticity' being of purely entropic origin. 
The deformations with q = 1 are, at least to linear order, translations of 
the cell with respect to the origin. As is briefly discussed in section 15.11 the 
expansion actually leads to the full probability density functional, valid in 
the limit n — > oo, of the zero-cell perimeter. 

The possibility of constructing an asymptotic series by expanding around 
the regular n-gon implies that when n — > oo, the shape of the n-sided cell 
tends with probability one to a circle (which we show to be of radius R c = 
(na/2n) 1 / a ). Miles [23] was the first to state this property for the Crofton 
cell and to validate it by heuristic methods. 

The approach to circularity in the limit of large sidedness n as studied 
here resembles, but is nevertheless distinct from, many mathematical results 
derived in the limit of large cell size, with the size being defined in various 
ways. These two limits correspond to ensembles that are fully disjoint. 

The conjecture originally formulated by Kendall in the early forties says 
that the Crofton cell becomes circular when its area goes to infinity. Calka 
and Schreiber [25] proved the approach to a circle when the radius of the 
largest possible inscribed disk tends to infinity. Hug, Reitzner and Schneider 
[TO] showed a generalized version of Kendall's conjecture in any dimension d 
when the fc-dimensional volume of the cell, in the sense of Hausdorff (2 < k < 
d), goes to infinity. Hug et Schneider [22] considered the zero-cell of a general 
class of d- dimensional tessellations which includes the two-dimensional family 
studied in this work, but also includes cells ('polytopes') resulting from more 
general direction dependent Poisson hyperplane processes. These authors 
prove the approach to a limit shape when the size of the cell, as measured in 
any of a variety of ways, tends to infinity. 

The present work, in summary, considers a special subclass of the systems 
introduced in reference [22] and distinguishes itself from that work in that 
it takes a different large-cell limit and leads to an explicit expansion for 
the sidedness probability. Our expansion in negative powers of n is not 
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mathematically rigorous; but it is of a kind that in statistical mechanics 
commonly leads to exact results. 



2 A one-parameter family of line tessellations 

We consider a set of lines defined by projections R that are distributed in 
the plane M 2 with density 

p(R) = \R a ~ 2 , R > 0, a > 1, (2.1) 

where A has the dimension of an inverse length to the power a and will be 
kept only as a check on the dimensionalities of our formulas. The expected 
number J\f(a, L) of projections in a disk of radius L centered at the origin is 
equal to 

f 1ir\T a 
Af(a,L)= dRp(R) = — — . (2.2) 
Jr<l « 

In order to have a well-defined problem in the infinite plane we first consider 

this disk occupied by N projections distributed independently according to 

the probability density 

P < R > = 2^(sf"' ° <R<L - (2 ' 3 > 

At some suitable point below we will take the limit L, N — » oo with N = 
Af(a,L). 



2.1 Sidedness probability p n (a) of the zero-cell 

Each projection R a , for a = 1, 2, . . . , N, has a line a associated with it. Then 
p n is the probability that n lines from among the iV contribute a segment to 
the perimeter of the zero-cell (that is, the cell enclosing the origin), and that 
the other N — n lines do not contribute segments. We will labeling these two 
groups of lines by a = 1, 2, . . . , n and a — n + 1, n + 2, . . . , N, respectively, 
and abbreviate R' n ^ = {Ri,...,R n }. Omitting explicit indication of the 
dependence of p n (a) on N and L, we can write 

p n {a) = ( N ) [ dRiP(Ri).../ dR n P(R n )x(R W ) 

\ n / JR!<L J R n <L 

xexp(-A*(R W )), (2.4) 
where x selects the R'"' that define Crofton cells: 

1 if Rl n l represents an n-sided cell enclosing the origin, 
in all other cases; 

(2.5) 



x(R [ ^ 
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Figure 2: Figure illustrating the definition of the support function h(&) given in 
section 12.1.11 

and exp( — ^4 a (R' n J)) is the probability that the other N — n lines do not 
intersect the perimeter of this cell. We will employ polar coordinates and 
write R a = (R a , $ a ). Upon using (12.31) for P and partially taking the limit 
N — > oo we can rewrite (12.4ft as 

Pn(a) = -r / d$i . . . d$„ / di?!^" 1 . . . di4#r 1 X(R W ) e-^( RW ). 

(2.6) 

In order to make progress we must now render A a and x explicit. 
2.1.1 Expression for A a 

Since the N — n lines are independent, exp(— A a ) is an (N — n)th power and 
we need to consider only a single line. We refer now to figure [2j For a half-line 
originating in the origin O and having an arbitrary angular direction $, let 
£($) be the tangent to the cell that intersects this half-line perpendicularly. 
Generically this tangent will have only a single vertex in common with the 
perimeter. The distance h(<&) from the origin to is called the 'support 
function' [26] of $. The no-intersection condition can now be restated as the 
condition that there be no projection R a (a = n + 1, . . . , N) in the sector of 
width A<3> within a distance h(<&) of the origin. The probability u;($)A$ for 
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a given R a to be in this sector is 



w($)A$ = / d$'/ dRRP(R) 



A$ //i($)y 
~2Tv Z J 



0(A$ 2 ), (2.7) 



where we used ( 12. 3ft . The probability exp(— ^4 a ) for the line through R a not 
to intersect the cell is equal to one minus the integral on all <3> of expression 
( 12.71) . Upon raising this result to the (N — n)th power and taking the limit 
N — > oo we obtain the desired weight exp(— A a ) with A a given by 



A ^ 



a 



A a (R w ) = - / d$ //*($). (2.8) 
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Equations (12.61) and (12. 8ft can also be obtained from Slivnyak's formula by 
the method of reference P3]. Formula ( 12. 6ft was originally included in [27]. 
Whenever (12.81) appears below, we will set A = 1, which entails no loss of 
generality. 

It is easy to find the explicit expression for the support function h(<&), 
which follows directly from its above definition. By a permutation of the 
labels of the polar angles and a rotation of the coordinate system one can 
always arrange things such that = $o < < • • • < ®n-i < 27r, where 
by convention $ n = <£> . We define the S m = (S m , \l/ m ) as the consecutive 
vertices of the cell in polar coordinate representation; in particular, S m is the 
intersection of the (m — l)-th with the m-th line (see figure [3]). After that a 
little algebra yields the two alternative expressions 

= [sin($ m -$ m _ 1 )] _1 [ J R m _ 1 sin($ m -$) + J R m sin($-$ m _ 1 )] 

= S m cos(^ m - $), $ m _i < $ < $ m , m = l,...,n (2.9) 

When one substitutes (I2.9P in (12.81) . a $ integral appears which can be carried 
out in closed form only for integer a. However, the large-n expansion to be 
introduced in the course of our development will apply for arbitrary a. 



2.1.2 Angular variables 

Other angles essential to our study may be defined in terms of the $ m and 
ty m and are shown in figure O First of all, 

£m = &m — $m-l , 

Vm = ^m+l-^rn, 771 = 1, 2, . . . , 71, (2.10) 

with the conventions $ n = $o + 27r and \P n = ^/q + 2tt. The £ m are the angles 
between two consecutive projection vectors and the r/ m those between two 
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consecutive vertex vectors; n-periodicity in their index m will be understood. 
For fixed sets of angles £ = {£ m } and r\ = {r] m } one may still jointly rotate 
the vertex vectors S m with respect to the projection vectors R m , as this 
modifies only the relative angles f3 m and 7 m (see figure [3]) between the two 
sets. We may select any one of these relative angles and call it 'the' angle of 
rotation, since it will determine all others. We will select (3\ for this purpose; 
when it is given, the remaining (3 m and 7 m can be expressed as 



m—1 



(3 m = Pi - ^(6 - Vl), m = 2, 



1=1 

m—1 



lm = -(3i + 2_^(€l-r}e)+€m, m = l,...,n, (2.11) 

i=i 

We now show that f3i cannot be arbitrary but is in fact determined itself by 
the two sets £ and rj. As is clear from figure [31 one can relate R m to R m -i 
by R m = (cos 7 m / cos/? m )i? m _i. Upon iterating n times starting from any of 
the Rm, periodicity imposes that we recover the initial value. Let us define 
Ghy 

e 2irG(tM = fT COS 7m 

J- J- cos/3 m 

m=l 

where the notation expresses that the 2n — 1 variables 71, 02, 73, . . . , fl n on 
the right hand side should be viewed as the functions (12.111) of the £ m and rj m 
and of the 'angle of rotation' f3\. Because of periodicity (3\ must then have 
the special value (3i = /?*(£, 77) that is the solution of 

G(£,r/;&)=0. (2.13) 

It was shown in reference [2U| that the solution of (I2.13P is unique. 

Finally we remark that whereas the £ m and r\ m must be positive, the f3 m 
and 7 m may have either sign. 

In terms of the angles defined above, equations fl2.8p - fl2.9l) can be made 
fully explicit in those cases where it is possible to carry out the integral on 
$. One such case occurs for a = 1. By Cauchy's integral formula we have 
A± = V, where V is the cell perimeter and has the two alternative expressions 

V = ^ Rm [ tan l£m + tan f £ m+ i] 

771=1 

n 

= Rm \ tan j m + tan /3 m +i] ■ (2.14) 

771=1 

For a = 2 we have that A2 is the area of the fundamental domain of the 
Voronoi cell, that is, of the union of n disks of radii S m centered at the vertices 
S m - Explicit expressions for A2 may be found in references [T6l 1281 [T7 1 fT8 l 120] . 
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■m—1 



Figure 3: Heavy line segments: the perimeter of the zero-cell in a random line 
tessellation of the plane. The R m are the projections of the origin onto the lines 
containing these perimeter segments. The S m are the cell vertices. The figure 
shows the angles £ m , rj m , (3 m , and 7 m defined by (|2.10p and (|2.1ip . 

2.1.3 Expression for the indicator function \ 

Having denned the angular variables we return now to the conditions imposed 
on the domain of integration in (12.61) by the indicator \. In terms of the angles 
(3 m and 7 m these conditions simplify greatly and take the explicit form [T7JHB] 

-f</#m,7m<f, /? m + 7m>0, 7 m + (3 m+1 > 0, (2.15) 

where m = 1,2, ... ,n and n-periodicity in m is understood. As could be 
expected, this condition depends only on the angles of the problem and not 
on the length scale. 

2.2 Transforming the expression for p n (a) 

Equation (I2.6P represents p n (a) as a 2n-fold integral which we will transform 
in successive steps to a form manageable in the large n limit. An important 
variable is the 'average radius' i? av defined by 



We now pass from the set of variables of integration {R m , &m} used in (12.61) 
to a new set of variables consisting of the radial scale i? av and the set of angles 
{Cm, Vm}- Although the (3 m and 7 m may be entirely expressed in terms of the 
£e and r)e, they constitute a useful set of auxiliary variables. 




(2.16) 



m=l 
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Convenient shorthand notation will be 

-27T />2vr 



\ = dfi & • • ■ d£„ / dr/i . . . d?7 n 5 ( ^ f m - 2tt) 5 ( ^ Vm ~ 2tt) 
Jo Jo m=1 m=1 

(2.17) 

and 

9 = J] 0(f - /3 m ) J] 0(f - 7 m), (2-18) 

m=l m=l 

where 9 is the Heaviside unit step function. After doing the algebra of the 
transformation of variables [TBI [20] we find that (12. 6p may be recast in the 
form 

Pn (a) = - [ G'&mfry^flfaT^-^e rWiC-'e-^, (2.19) 

in which G' denotes the derivative of 77; j3\) with respect to fti] we have 
abbreviated T m = sin£ m / cos 2 /3 m ; and the 2n ratios p m = R m /R w may be 
expressed in terms of the angles by means of the 2n relations 



p m = (cos7 m /cos/5 m )p m _i (m = 2,...,n), n 1 ^p m = l, (2.20) 



m=l 



of which the last one follows from (12.161) . 

Clearly A a can be written as R" v times a function of the angles. We will 

set 

Aa = — R^ + n-'V), (2.21) 
a 

where 2ira~ l R® v is the value that A a takes for a circular cell of radius _R av 
and in which V depends exclusively on the angles. The notation in (12.211) is 
meant to suggest that n~ x V is of order n" 1 as n — ► 00, but we will rely on 
that only later. Integrating (12.191) on R av converts it into 



. . (n - 1)1 f a \ n ~ l 
Pn{a) = - 



(srX 66- ^' a - 1 ' (2 ' 22) 



2-irn \ 2ix 
where 

n 

e" V(Q) = G\t V, [ II P™ T m&] (1 + n-HO"". (2.23) 

m=l 

The integrations that remain on the right hand side of (I2.22p bear only on 
the angles £ m and r] m , that is, on the shape of the cell without regard to its 
radial dimension. We still remark that the factors in the product on m 
in (12.231) compensate the factors £ m incorporated in the definition (I2.17P of 
. Their purpose is to ensure that V(a) remains finite when any of the £ m 
tends to zero. 
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We finally write (12~22D as 



Pn(<*)=P$ ) («)(6e- v ^) 



a > 1 



(2.24) 



where for any function X of the angular variables we define 



(X) 




(2.25) 



Straightforward calculation yields 




(n 



i V2tt 



) 



n-1 




1 



27m 




n-1 



a > 1, 



(2.26) 



which is equation (jl.2p of the introduction. Equations d2.24[) and (12.261) 
represent an important step forward with respect to the initial expressions 
( 12.61) and ( 12. 81) . However, the hard part of the problem remains, namely 
to determine the n dependence of (Qe~ Y ^) in equation (I2.24p . In sections 
l2T4H231 we will show that for a > 1 lim n _ >00 (6e- v ( Q )) = C(a), where C(a) is 
a finite numerical constant. The special case of the Crofton cell, a = 1, will 
not be covered by the arguments of those sections and must be dealt with 
separately; we will do so in sectional 

2.3 Scaling of the average radius i? av 

The inner integrand of (I2.19p . when combined with (12.211) . shows that the 
average radius i? av has a probability distribution 



which depends on the cell shape through V. The notation is meant to suggest 
that n~ 1 V is negligible when n — > oo, and this will be confirmed below. 
Elementary analysis then shows that for n — > oo the distribution P av (i? av ) 
has a peak at R c whose width is a c , where 



Pav(^av) = CSt X R ^n-X ^MR^l+n 



a > 1, 



(2.27) 




(2.28) 



Hence the typical deviation of R av from R c is a factor n2 smaller than R c 
itself, for all a > 1. 
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2.4 Large-n expansion of Y(a) 

The regular n-sided polygon is a point of high symmetry in phase space where 
we expect the integrand of in f 1 2 . 2 2 j) to be stationary. This point has 
£,m = Vm = 2nn~ 1 and (3 m = ^ m = im~ l for all m = 1, 2, . . . , n. Coordinates 
that describe the deviations from this symmetric state are defined by 

$£m = £m - 27rn _1 , 5rj m = rj m - 2im~ 1 (2.29) 

and their suitably scaled Fourier transforms 

in In 

X q = ^Y. e 2niqm/n ^m , % = ^ e 2 ^ m/n S Vm , (2.30) 

m=l m=l 

where q = ±1, ±2, .... It will appear that the integrations required for the 
calculation of the average in (12.241) can be done explicitly once expressed in 
terms of the variables of integration X q and Y q . 

2.4.1 Scaling with n 

We will treat V perturbatively in inverse powers of n. An initial hypothesis 
on the smallness of the angles with n is suggested by the scaling that prevails 
when in (12.221) we set V = 0. That yields 

££m, Hn ~ n' 1 , (2.31) 

where the symbol ~ indicates the scaling with n in the large-n limit. Using 
relations (I2.1ip and (I2.20p between the angles as well as the fact that the £ m 
and rj m are to leading order independent, we find from (12.311) that furthermore 

/3 m ,7 m ~n~5, r m = p m - 1 ~ n~2; (2.32) 

and using (I2T3HD we find from fl2T2^|) that 

X q ~ n°, Y q ~ n°. (2.33) 

If and when needed, the /3 m , j m , and r m and their Fourier transforms can be 
expressed in terms of the X q and Y q . 

The initial hypothesis now holds that the above scalings remain valid in 
the presence of the nonzero V defined by (I2.23p . We will consider the validity 
of this hypothesis confirmed when at the end of this section we will find that 
in the limit n — >• oo it produces a finite leading order result for (exp(— V)). 

We end these considerations with two remarks. First, the present scalings 
are identical to those encountered in the study of the large n-sided Voronoi 
cell. How to handle them technically, in particular when sums of ~n terms 
appear, is nontrivial and has been discussed in detail in reference [H], sections 
5 and 6. Secondly, we note that the smallness of the r m with n implies that 
all Rm are close to the average R av , which in turn together with the smallness 
of the £ m implies that for n — >• oo the perimeter approaches a circle. 
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2.4.2 The expansion 

We defined V recursively by (12T23]) . f[2~2ll . and (JZB]) in terms of V, A a , 
and /i( < &). Since h(<&) defined by (12. 9p is piecewise analytic on n successive 
angular intervals, when substituted in (12. 8p it gives rise to a sum of n terms. 
In the mth term, that is, for $ m _i < $ < $ m , we set S m = R m /cos~f m and 
pass to the variable of integration 9 = $ — \P m . Then (12.81) becomes 



A a = -Y(-^) / d#(cosi 



™— 1 1 9 fm ~r • • • J -0m 



Here the dots indicate terms of higher order in the angles. We may expand the 
ath powers in (12.341) and do the 9 integral. Using the sum rule Y^m=i (Pm + 
7 m ) = 2n we find that the leading order result for A a is (27r/a)i?° v . Com- 
parison to (I2.2ip shows that the higher order terms in (I2.34p determine V. 
Pursuing the expansion to higher orders we find 

2tt ™ - 

— V = a ^ r m (7 m + /3 m ) + ~a(a r^(7 m + (3 m ) 

^ i 1 

m=l m=l 

+ \a{2rf m + 3 7 ^ m -(31) + .... (2.35) 

Equation (12.351) for V may be substituted in (12.231) . The other factors in 
(I2.23P may be expanded similarly, and together this leads to the expansion 
of V. The order in n of each term in the expansion may be estimated in 
the way outlined in section 12.4.11 It then appears that Y(a) allows for an 
expansion in powers of n~2 ; 

V(a) = Vi(a) +0(tT3), (2.36) 

of which the leading term Vi(a) is a quadratic form in the angles that is of 
order n°. Upon expressing all variables in terms of the X q and Y q one obtains 

Vi(a) = J2(X g ,Y q )-V q -(X. q ,Y. q ) T , (2.37) 

where the superscript T indicates transposition and where V g is the sym- 
metric matrix 

V q =\ " 2 2.38 

\-A q + \B q A q -B q ) 

with 

qz 2g 4 q z 

Equations (I2.36I) - (I2.39I) complete the large-n expansion of V(a). They are 
valid for all a > 1. 
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2.5 Large-n expansion of p n (a) 

The large-n expansion of p n (a) is based on the one of V given above. Upon 
substituting fl2.36p - fl2.37p in (12.241) and using that 6 may be replaced with 
unity up to corrections that vanish exponentially for large n, we get 

p n (a)~p£°)(eO<e~ Vl(a) ) (2.40) 

with Vi given by (12T37D . 

Although (12.371) is a quadratic form, the X q and Y q are not Gaussian 
distributed. Nevertheless, it was shown in reference [H] that to leading 
order in n~ 1 / 2 the average (. . .) in (I2.40p may be carried out as though the 
X q and Yg were Gaussian, with a probability distribution 

M G exp (- Y,{K%)-V- X -{X_q,Y_ q Y), (2.41) 

q^O 

where Mg is the appropriate normalization constant and where we introduced 
the 2x2 diagonal matrix E = diag{l, 2}. This fact is partly confirmed by the 
result of convergence in distribution and asymptotic mutual independence of 
the Fourier coefficients of i.i.d. random variables [2^1 [3U]. The Gaussian 
integrals are easily done. Still letting 1 represent the 2x2 unit matrix one 
finds that 

oo 

p n (a)~C(a)p£\a), C(a) = J] A? (a), a > 1, (2.42) 

where we abbreviated for all a > 1 

A q (a) = det(l + V g E) 

= 1 + 3A q - 2B q - \B 2 q 

= 1 - (3 - a)q~ 2 + a 2 q~ A . (2.43) 

Equations (12.421) and (I2.43P are the main result of this section. It is easily 
verified that A q (a) > for all a > 1 and q = 1, 2, . . . with the only exception 
that Ai(l) = 0. This is the reason why the substitution that led to (12.401) 
must be restricted to a > 1; for the Crofton point a = 1 the resulting 
Gaussian integration would diverge. Hence the method of this section leaves 
the Crofton problem unsolved. In the following section we will develop a 
modified method suitable for the Crofton cell. 

3 The Crofton cell: a = 1 
3.1 Centering condition 

For a = 1 the expansion of (Qe~ Y ^) around a regular n-gon leads to a di- 
verging integral and hence cannot be used for finding p n (l). This divergence, 
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due to the vanishing of A 9 (l) with q = 1, may be tracked down to the ori- 
gin being located with the same probability in any point of the cell; hence 
the problem has a 'zero mode' due to translational invariance. Indeed, for 
q = ±1 elastic deformations of a circle are actually translations. 

This analysis also points the way out of the difficulty In the integrals 
one should count a Crofton cell only when, under a suitable definition, it is 
'centered' around the origin, and then give it a weight proportional to its 
area in order to simultaneously account for all noncentered cells that may be 
obtained from it by translation. The 'centering condition' that we adopt is 
that the q — ±1 Fourier components of R m vanish, that is, 

1 - 

R x = - e 2nim/n R m = 0. (3.1) 

m=l 

This is a constraint only on the vector lengths R m , irrespective of the angles 

The following remark establishes one more link with the mathemati- 
cal literature. There, a function *4. 2 (R' n ^ r ) appears defined as A a (R} n ^) 
with a = 2 in (I2.8P except that the support function /?.($) with respect 
to the origin O is replaced by the support function /i(<3?; r) with respect 
to r. The function *4. 2 (R^; r ) is commonly used to calculate the proba- 
bility that a convex set containing the origin translated by — r is included 
in the typical Poisson-Voronoi cell. As a function of r, «4 2 (R j r ) is a 
convex function such that its gradient V r .4 2 (R^;r) is, to, up to a multi- 
plicative constant, equal to ( /i($; r) cos $6?$, J Q 27r /i(<3>; r) sin $<i$) . The 
sum appearing in (13. 1 p is the discrete version of this gradient. Hence (13.11) 
amounts to adopting as the 'center' of the cell the point for which the fun- 
damental domain has minimum area; we remark that this point is given by 
(tt- 1 J q 2 " h{$) cos $ d$, 7T- 1 J Q 2 " h{<!>) sin $ d$) . 

In the presence of constraint (13.11) the necessary mathematics requires a 
few preliminaries. Let t = (t x ,t y ) be a translation applied to a Crofton cell. 
Under t the vertices S m are translated by just that amount. However, the 
R m , being projections of the origin onto lines that under the translation stay 
parallel to themselves, transform differently. If t keeps the origin inside the 
cell, the R m transform into Rm,t = (Rm,t, &m,t) given by 

cos $ m + t y sin $ 

$m,t = $m- (3.2) 

We will write = {Ri,t, • • • , R n ,t}- It will be convenient to extend the 
definition (I3.2p of R^ to translations t that take the origin out of the cell. 
In that case, however, R^ is 'nonphysical' (it does not represent a Crofton 
cell any more; one or more of its R m ,t are negative) and, by virtue of (12. 5p . 
the indicator x(Rj^) then vanishes. 
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We now show that if Ri n l defines a Crofton cell, then there exists a unique 
'centering translation' t*(R^), that is, one for which (13.11) is satisfied: 

Ri,r = 0. (3.3) 

The proof goes by explicit construction: using (13. ip and (13. 2p in (13.31) leads 
to a linear system of equations for t* whose solution is 

r(RN) = _I^^^, (3.4) 



where for any complex variable z we set z = (z c + iz s )/y2 and where if is 
the matrix of elements 

(fu = cos(27rm/n) cos $ m , <y2 12 = cos(27rm/n) sin $ m , 

m m 

(f2i = sin(27rm/ri) cos $ m , V?22 = ^""^ sin(27rm/n) sin $ m . (3.5) 

m m. 

Two useful relations, 

t*(R[ nl ) =t*(R M ) -t, (3.6) 

and 

/ dt5(t-t*(RM)) x ( R N) = l. (3.7) 

will serve below. 



3.2 Modified starting point for a = 1 

After these preliminaries we return to the general expression (12.61) for the 
sidedness probability p n (&)- For the case of the Crofton cell we insert (13.71) 
and rewrite p„(l) as 

\n /"27r p poo 

Pn(l) = -T d$!...d$ n / dt / d^.-.d^ 

™ Jo ii 2 Jo 

x *(t - t*(RN)) x(R N) x(R W) e-'CaM), (3 _ g) 

where we used that «4i = P and have indicated explicitly the dependence on 
RM of all quantities involved. 

Under the t integral in (13. 8p we now transform from the R m to new radial 
variables of integration R' m = R m ,t, with the R m ,t given by (13.21) . We may 
express R^ in terms of R'M by means of RW = R'J™ 1 . The Jacobian of 
this transformation is unity but the limits of integration should be treated 
with some care. Equation (13.21) implies for R' m the domain of integration 
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t x cos $ m + ty sin $ m < i?^ < oo. If the lower limit of integration is negative, 
we may replace it with zero because x(Rt^) = on the interval discarded; 
and if the lower integration limit is positive, we may also replace it with zero, 
because then x(Ri n ') = on the interval added to the domain of integration. 
Hence, still using (13 .6p . we get 

\n p2ir p poo 

Vn{l) = -r d$!...d$ n / dt / dR[...dR' n 
n\ Jo Jm.2 Jo 

x 5(t*(R'W)) X (R'N) x(R 'N) e -v(n'™)_ (3.9) 

Since the cell perimeter is invariant under translation, we have P(R_™ ) = 
V(R!^). The only t dependence left in the integrand of (I3.9P is then the 
one in x(R'J^). The integral on t of this quantity has a nonzero contribution 
only in a domain of the same size and shape as the cell itself and therefore 
produces the cell area A(R'^). Hence, still suppressing the primes, we get 
from ([33D 

\n /■27T poo 
Pn (l) = — d^... d$ n dRL.-dRn 

™ Jo Jo 

x 5(t*(RW)) x(R [n] ) A(RN) e ^ RM \ (3.10) 

Equation ( 13.101) constitutes a modified starting point for the calculation of 
the sidedness probability p n (l) of the Crofton cell. It differs from the original 
expression (12. 6ft by the insertion of a delta function and of the area factor 
A(RN). 

3.3 Sidedness probability p n (l) of the Crofton cell 

Having rewritten the definition of the sidedness probability p n (l) as (I3.10p . 
we are now in a position to start its explicit evaluation. The cell area A is 
given by 

n 

A = R^A^g, A ang = |^ (tan 7 m + tan/3 m+ i), (3.11) 

m=l 

where A ang is an expression only in terms of the angular variables. We may 
similarly isolate the radial part of the centering translation t* found in (13.41) 
by writing 

t* = RlXn g i (3-12) 

where t* ng depends only on the angles. 

When we substitute (13.111) and (13 . 1 2[) in (I3.10p . the factors R av stemming 
from 5(t*) and from the area A cancel. We set V = 2irR av (l + n~ 1 V), which 
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is the special case a = 1 of (12.211) ; here 27iR av is the perimeter of a circular 
cell of radius i? av . We may in the same way as before pass to the variables 
{£,m,Vm} and a single length scale _R av . Carrying out the _R av integration we 
get, setting A = 1 as before, 



-v 

ang X e 



2 i*z7v poo 
Pn(l) = -T / d$!...d$ n / d^L.-di^t:^)^ 

n! Jo Jo 

If r n i f 00 

= -/ CCe.i?;/^)- 1 Ip^" 1 e / diCH 



n-l 



av - 1 i-av e 



(n-l)! 



5(t: ng )A ang e- v ^, (3.13) 



n(27r) n 

which is an alternative for (I2.22p when a = 1. We may rewrite (13.131) as 

Pn(l) =Pi 0) (l) (e5(t; ng ) A ang e- V «), (3.14) 

where and V(l) are given by the same equations as before, namely 

( I2.26j) and (I2.23f) with a — 1, and the angular brackets are defined by (12. 25ft . 
When a — 1, expression (13. 14f) is an alternative to (12 .24j) . Both are exact 
representations of the initial expression for p n (l), but the (13. 14j) is required 
if the large-n expansion is to succeed. That will be the subject of section |2~41 

3.4 Large-n expansion of p n (l) 
3.4.1 Preliminaries 

We will perform the large-n expansion of the average appearing in (I3.14p . 
The expansion of V(l) is not different from the general case and the result 
follows directly from (I2.36p - (l2.39p by setting a = 1. As before, the indicator 
G will give corrections that are exponentially small with n; we may therefore 
replace it with unity. The insertion v4 ang in (13.141) . given explicitly in (13. lip , 
is expanded as 

n 

^ang = [lm + Pm+1 + 0(7™X + l)] 

m=l 

= 7T + C>(rH), (3.15) 

which, to leading order, is the area of a disk of unit radius. We now investi- 
gate the factor <5(t* ng ) for large n. In view of (13 . 12[) and (13 ,4p this requires 
the asymptotic evaluation of the matrix ip defined in ( I3.5p . In the large-n 
limit the angle differences ^ = $^ — §e-i are independent random variables 
of average 2-ir/n. Hence after setting $ = and summing ^ on I from 1 to 
m we have $ m = 2irm/n + OirT^\ By doing the sums in (13. 5p one then 
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finds that <*p is equal to \n times the 2x2 unit matrix, up to corrections 
that vanish for n — > oo. Setting R\ = R av pi we have 

where the symbol ~ indicates validity for n — > oo. In that limit we can 
express pi in terms of the X q and Y q by means of the relation [IB] 

ft^n^Xi-Yi). (3.17) 

Using fl3TT5|) - fl3TT71) in (jUIJ) we find that 

p n (l) ~ pf(l) x §7m - ^ C )(J(X 1 8 - n s ) e- Vl «), (3.18) 

where Vi(l) is given by (I2.37P with a = 1 and (. . .) is the average with 
respect to the Gaussian weight (12.411) . 



3.4.2 The Gaussian integrations 

In (I3.18P the integrations on the Fourier variables X q and Y q with q ^ ±1 
may be carried out in the same way as before and lead to 

oo 

Pn (i) ~ ^(i)xi TO [i(i-f 2 r 2 

q=2 

x (6(Xl - Y*)8{Xl - YD e -*i*-i+*i*'-i) j (3.19) 

where the terms in the exponential represent the q = ±1 contribution to 
Vi(l). The average (...) in (13.1911 reads explicitly 

(...) = ^2 / dXidX^dYidfli . . . e~ 2jt ^-^-\ (3.20) 

as follows from (12.411) when restricted to its q = ±1 Fourier components. 
When (I3.20p and (13.191) are combined, the terms Y{Y_\ in the exponents 
cancel and the average in ( 13. 19ft turns out to be equal to l/(37r). Without 
the delta function insertions in the angular brackets, these final integrations 
on Kj-i would have led to a divergence. The fact that they now remain finite 
confirms the validity of our approach. 

After combining everything and still using that n^l 2 (l — q~ 2 ) = |, we 
obtain the main result of this section: the sidedness probability p n (l) of the 
Crofton cell is given by 

Pn(l)^|npi 0) (l), n^oo, (3.21) 

with Pn (1) given by (11.21) . This is what was announced in section [l~4l of the 
introduction. 
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4 The typical cell in a Poisson line tessella- 
tion 



In this section we continue our study of the case a = 1, that is, the Poisson 
line tessellation. But whereas in section[3]we studied the Crofton cell (or zero- 
cell), we will now consider the typical cell. The preceding results give rise to 
a corollary concerning the typical cell that can be derived with little effort. 
Let p typ (l) denote the sidedness probability of the typical cell in the Poisson 
line tessellation. As was shown in detail by Calka p3], the expression for 
p typ (l) differs from equation (I3.8p for p n (l) by the insertion in the integrand 
of the latter of an extra factor (A) typ / A(R^) , where (A) typ = 1/vrA 2 is the 
average cell area [2] in the Poisson line tessellation. 

The extra l/A(R' n J) cancels the A(Ri n l) present in (I3.1(jp and hence 

Pn yP (l) = ^—r d*i ■ . .d$ n / dR 1 . . .<LRn<f(t*(RM)) X (R W ) e^ (RW) . 

TTTl! J Jo 

(4.1) 

In the limit n — > oo this is easily evaluated by the methods of sections 13.31 
and 13.41 The integral on i? av now has an extra factor in its integrand 
and as a consequence we get instead of (13.141) the expression 

4.77- 

PTW = 7 7y 7y\ (©^tLg) (1 + n- 1 y)- 2 e- v «), (4.2) 

(n — l)(n — 2) B 

where we have again set A = 1. Equation (14.21) is still a fully exact expression 
for the sidedness of the typical cell. ^From here on, again, we have to resort 
to a large-n expansion. Since (l + ri" 1 ^)" 2 = l + C^n" 1 ), the average in (14.21) 
is to leading order in n identical to the one in (13.141) except for the absence 
of the insertion A ang , which to leading order is equal to n. Upon combining 
these considerations we get from (14.21) 

p**(l) ~ §71-^(1), n^oo. (4.3) 



This is our final result, announced in section [T~4"l for the asymptotic sidedness 
probability of the typical cell. Comparison of (13.211) and (14. 3 p shows that the 
zero-cell and the typical cell sidedness probabilities have the ratio 

where 7ri? 2 = n 2 /47r [see (12.281) ] is the average area of an n-sided Crofton 
cell and (A) typ = is the average typical-cell area. Equation (I4.4p is 
heuristically obvious. 
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5 Properties beyond p 



5.1 General 

In this section we return to a general value of the parameter a. The analysis 
of sections [2] and [3] focused on finding the sidedness probability p n . It started 
from a weight functional for the zero-cell expressed initially [in ( 12. 4ft ] as 
X exp(— A a ). In the course of the analysis we were led to transform this func- 
tional into exp(— Vi) where Vi is quadratic in the Fourier variables X q and 
Y q , which in turn are Gaussian distributed according to (I2.4ip . The transfor- 
mation is valid in the limit n — > oo and the resulting weight exp(— Vi) may 
be used to calculate arbitrary averages of functionals of the cell perimeter, 
provided that sums and products on the wavenumber q converge sufficiently 
rapidly. We will limit ourselves below to discussing only one of these cell 
averages. 



5.2 Fluctuations of R m around R av 

Important information on the cell shape is contained in the root mean square 
deviation <tr of the vectors R m from their average i? av . Straightforward 
calculation [T8] yields for this quantity the expression 

n 

4 EE n^YsdRm-R^) 2 ) 

m=l 

oo 

~ 2n- l RlY,q-\\X q -Y q \ 2 ), n^oc. (5.1) 

9=1 

It only requires a Gaussian integration to find that for n — > oo 

(\X q -Y q \ 2 )~lA-\a), g = 2,3,..., a > 1. (5.2) 



A special case is 



i^-nn-^ 1 ; (5.3) 

U, a = 1, 



where the discontinuity at a = 1 is due to the centering condition that we 
applied in that case. Substitution in ( 15. ip of ( 15. 2ft and ( 15. 3p . as well as of the 
explicit expression ( 12. 28ft for R c , shows that 



j__ i 

cfa) n a 2 , a > 1, 



c(l) n 2 , a = 1, 
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where the coefficients have the properties 

1 

— T, Oi— ► 1, 

2n(a - 1)2 

00 

= ^E(^Tj^ = ^( 1 - gr) = 0.010256... (5.5) 

The divergence of c(a) as a — > 1 is due to the divergence of the q = 1 term 
in the sum in (15.11) ; it demonstrates once more the tendency of the center of 
the cell to become delocalized from the origin in that limit. For a = 1 we 
obtain discontinuously the value of c(l) due to our centering of the cell. 

6 Conclusion and final remarks 

We have determined the asymptotic large-n expansion of the sidedness prob- 
ability p n (ct) of the zero-cell in a family of tessellations dependent on a pa- 
rameter a. Special cases are the typical Poisson-Voronoi cell (a = 2), and 
the Crofton cell (a = 1). The latter turned out to be a singular point in the 
parameter range and required an extension of previously known methods. 

We end with a few remarks about further research to which this work 
opens the door. The first remark concerns the parameter range < a < 1. 
In this range the density p(R) and the probability distribution P(R) defined 
in (12. 3p are still integrable in the origin. However, Ai(a) is negative, which 
indicates a tendency for the origin to be away from the cell center defined by 
equation (13. 11) . Exploring the range < a < 1 further would definitely be of 
interest; it is no longer certain that in that case, too, the zero-cell will tend 
to a circle in the large-n limit. 

Secondly, the analytic method of this work points the way to a new Monte 
Carlo simulation algorithm for the Crofton cell, similar to the one that was 
used [20] successfully for the Voronoi cell, and capable of producing accurate 
p n values for all finite n. 

Thirdly, one may consider sidedness correlations between distinct cells, 
whether neighboring or separated by a larger distance. In the Voronoi case 
(a = 2), the correlation between adjacent cells is commonly denoted m n and 
defined as the average sidedness of a cell given that its neighbor is n-sided. 
To our knowledge, no one has yet attempted to determine the analog of this 
quantity for the Poisson line tessellation. 

We will leave these matters for the future. 



c(a) 
c(l) 
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